Wave functions of linear systems 



o 
o . 

CN ■ TOMASZ SOWINSKI 

^ . Center for Theoretical Physics 

of the Polish Academy of Sciences 
I/-) Al. Lotnikow 32/46, Warsaw, Poland 

tomsow@cft.edu.pl 

> 
O 

m 

Complete analysis of quantum wave functions of linear systems in an 
OO \ arbitrary number of dimensions is given. It is shown how one can con- 

■ struct a complete set of stationary quantum states of an arbitrary linear 

' system from purely classical arguments. This construction is possible be- 

cause for linear systems classical dynamics carries the whole information 
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i <~"| about quantum dynamics. 
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1. Introduction 



Connection between classical and quantum description of physical sys- 
tem manifests itself in very amazing and nontrivial way. We understand 
how to describe any system in classical and quantum language and also 
we believe that quantum description should smoothly transform to classical 
one when we neglect quantum correction^]. But in general we still do not 
understand how one can make this transformation in practice and identity 
the classical motion in quantum one. The first description how one should 
understand this correspondence was given by Ehrenfest in his famous theo- 
rem [I]. It states that the expectation values of physical observables evolve 
in time almost as classical quantities (up to quantum corrections). In the 
case of linear systems they evolve exactly in the same way. This is evident 
in the Heisenberg picture of quantum evolution. Since in the case of a lin- 
ear system classical Hamiltonian equations of motion are linear they have 
exactly the same structure as quantum Heisenberg equations of evolution 
of quantum operators. Therefore, in this case the classical dynamics carries 
the full information about quantum dynamics. 



1 It is often said that we go with h to zero. 



(1) 



2 



WAVE 'FUNCTIONS PRINTED ON FEBRUARY 1, 2008 



The question which we answer in this paper is how this exact connection 
between classical and quantum dynamics in arbitrary number of dimmen- 
sions is realized in the Schrodinger picture where the whole dynamics is 
contained in the evolution of the wave function of the system. Full dis- 
cussion of this correspondence in one dimensional case was given before in 

The main observation presented here is a multidimensional generaliza- 
tion of our idea proposed earlier in [3] for the three dimensional case. 

2. Class of linear systems 

For each N-dimensional linear system there always exist such canonical 
variables £ and tt that the system is described by the Hamiltonian of the 
following form 

I ~ /s TYY j 

H=—ir'F-ir + S'Q'n+-S-U'Z + mf(t)-S + -h(tyir, (1) 
2m 2 m 

where the matrices F and U are symmetric and the vectors f(t) and h(t) 
are given functions of time. If iv is chosen as a kinetic momentum then 
the matrix F should be positive definite (kinetic energy should increase 
with momentum of the particle). Therefore there exists a matrix O which 
diagonalizes the matrix F to identity 

T -F-6 = i. (2) 

Let us define the following matrices 

W = O^QO, (3a) 

S = ±(w + W T ), (3b) 

n = ±(w-w T ), (3c) 

and now let us make the following canonical transformation to the new 
canonical variables r and p 

r = 6 T -t (4a) 
p = d- 1 -Tv + mS-d' T -€+ f6 T ) _1 -/i(t). (4b) 



One can show that it is indeed a canonical transformation because the Pois- 
son bracket structure is preserved. This transformation is in fact the com- 
position of three simpler canonical transformations: the diagonalization of 
the kinetic energy term 

r' = O t • £, p' = Q- 1 ■ 7v, (5a) 
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the shift in the momentum space 

r» = r >, p" = p'+(0 T ) -1 -M^ 
and the stretching transformation of the phase space of the system 

r = r" , p = p" + mS • r" . 
In our new variables, the Hamiltonian of the system has the form 



p m 1 

H = h r il p H r-V-r — mg(t)-r 

2m 2 ' 2m 



where 



V 
9(t) 



6- 1 -u-(6 T 



s 2 



J 



■hit) 



1 -w-(d T ) l -h(t)-6-f{t). 



m 



(5b) 
(5c) 

(6) 

(7a) 
(7b) 



The last term in the Hamiltonian ([6]) depends only on time. Therefore, 
without loss of generality we can omit it. 



3. Classical dynamics 

Classical equations of motion following from the Hamiltonian © have 
the form 

^ = P-Qr, (8a) 

dt m y ' 



dp 

— = —mV-r — Q-p — mg(t). 



(8b) 



Of course the best way to solve these equations is to solve first a homoge- 
neous equation with g(t) = 0. In that case any solution of equations (jSJ) 
can be expressed as a sum of the eigenmodes 



r(t) 
P(«) 



(9) 



where the coefficients are determined by initial conditions. The char- 
acteristic frequencies and the amplitudes (Rk,Pk) obey the following 
matrix equation 
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In the Appendix |A] we show that for such a matrix, the eigenvalues appear 
in pairs. It means that for an N dimensional system we always have the 
following frequencies 



11 



This means that amplitudes of the modes tUi and — tUi are related by complex 
conjugation. Therefore, any physical solution of Eqs. © can be expressed 
in the following way 



Tit) 
P(t) 



E 

i=l 



R, 



p* 



-iujit 



(12) 



In the presence of the external field the only difference is that the parameters 
Aj are time dependent and the general solution has a form 



r(t) 
P(t) 



E 

i=i 



Ai(t) 



iuJit 



m 



(13) 



Time evolution of the coefficients Aj(i) is determined by the time dependence 
of the vector g(t) and is given by the ordinary differential equation 



d\j{t) 
dt 



-mgi(t), 



(14) 



where gi(t) are the coefficients in the following expansion of the vector g(t) 
into the mode amplitudes 





g(t) 



E 



9i(t) 



+9i(ty 



P 



(15) 



A full classical analysis of linear systems in three dimensions was given 
before [21 [1] . One can easily generalize this analysis to an arbitrary number 
of dimensions. 



4. Quantum dynamics of wave packets 

In the position representation of a quantum system, the wave function 
fy(r,t) obeys the following Schrodinger equation 

ihd t ^(r,t) = (-—V 2 + ^r-A-V + — r ■ V-r + mg(t) -r \ *(r,t). (16) 
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4-1- Dynamics of Gaussian wave packets 

In the first step of our analysis, let us consider the dynamics of a state 
of a quantum-mechanical system described by a Gaussian wave function of 
the form 



*(r,i) =N(t)e 



i<t>{t)/h 



cxp 



m 
2h 



(r-R(t)).K(t)-(r 



R(t)) + -r-P(t) 



(17) 

where the matrix K{t) is of course symmetric and its real part is posi- 
tive. This wave function should also be normalized. Since the Schrddinger 
equation (|16p preserves normalization in time, it is enough to assume the 
normalization condition at the initial moment t = 



1 = f d N r|^(r, 



0)| 2 = N(0f 



2irh 



mDet(UK(0)) 



N/2 



(18) 



Parameters R and P have a direct interpretation as the position and mo- 
mentum of the center of the wave function and at the same time they are the 
expectation values of the quantum operators of the position and momentum. 

One can show that the Schrddinger equation (fT6|) is equivalent to the 
following equations for the parameters of the wave function 



dK(t) 

dt 
dR(t) 

dt 
dP(t) 

dt 
dN(t) 

dt 

m) 

dt 



-ik{tf + iv- n,K(t) 
p(t) 



in 



n-R(t), 



-mV-R(t) - n-P(t) - mg(t), 



-Tr(3?K(t)), 



-Tr(KK(t)) - -±L + -R(t)-V-R(t). 



(19a) 
(19b) 
(19c) 
(19d) 
(19e) 



Comparing the equations (|19bj) and (|19c|) with the classical equations 
of motion ([8]), one can see that the dynamics of the center of the wave 
packet is the same as the dynamics of a classical particle. It is a well known 
realization of the Ehrenfest theorem which is, as we have said before, exactly 
satisfied for linear systems. 

It is worth while to notice that the dynamics of the center of our wave 
packet is totally separated from the dynamics of the internal motion of 
the packet. This is a general property of quantum dynamics in external 
harmonic potentials [Tj- In addition, the evolution of the shape of our 
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wave packet (described by a complex matrix K) is not influenced by the 
external field git). The only place where the external field appears in our 
equations is the dynamics of the center of mass. Therefore, to solve the 
general problem of the motion of the Gaussian wave packet, one can solve 
at first a simpler problem of quantum dynamics without an external force 
and without motion of the center of the packet. Then one just needs to 
move this solution along a chosen classical trajectory. 



4-2. Time evolution of the shape of the wave packet 

The equation ()19a|) describes the dynamics of the shape of the Gaussian 
state and it has the form of a well known matrix Riccati equation. Math- 
ematical theory of these equations is well developed [10] . In particular, 
the matrix Riccati equations find numerous applications in control theory 
[6] . It is known that these nonlinear equations can be replaced by the linear 
ones. Following the standard procedure [5], we shall search for solutions of 
Eq. (19a) in the form 

K(t) = --N(t)-D~ 1 (t). (20) 
m 

This matrix equation is satisfied (as is shown in the Appendix [X]) when the 
matrices N and D obey the following linear equations 

dN ~ - - 

= -mV-D - n-N, (21a) 

at 

— = —N — QD. (21b) 
dt m 

The linearization of the Riccati equation, in addition to being an effective 
mathematical tool, has also conceptual advantages. Namely, it leads to a 
direct relationship between classical and quantum theory. Comparing Eqs. 
(|2ip with Eqs. ([8]), one can see that the columns of the matrices N and 
D satisfy the same equations as the classical position and momentum vec- 
tors, respectively. Therefore, from the knowledge of the classical motion, 
one may determine the evolution of Gaussian wave function. It is a desir- 
able manifestation of an exact connection between classical and quantum 
mechanics in the language of wave functions. 



4.2.1. Example in one dimension 

Let us show now how this formalism works for the one dimensional 
standard harmonic oscillator [3 [8]. In this case our system is described by 
the Hamiltonian 

2 2 

p rntu 9 . . 

H = h + — *■ (22 > 
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In this case, the equations of the classical motion have a simple, scalar form 

. &H p . dH 2 
x = tt- = — , P = —-ft- = —mco x. (26) 
op m ox 

The solution of these equations satisfying the initial conditions x(0) = xq, 
p(0) = po is the following 

x(t) = xq cos(wt) + — ^- sin(a;t), (24a) 
muj 

p(t) = po cos(cot) — mujxo sin(o»t). (24b) 

In the quantum mechanical case, the dynamics is described by the Schrodinger 
equation 

/ h 2 d 2 \ 
ihd&fat) = - — — I +muj 2 x 2 V(x,t). (25) 



2m dx 2 

Let us assume that the wave function of the system has a Gaussian form 

where the parameter a(t) describes the shape of our Gaussian packet. Ac- 
cording to the equation (|19ep , the phase of our wave function is determined 
by the evolution the of a(t) 



1 



4>{t) = - j J o »a(r)dT. (27) 

If we know that at the beginning the shape of our packet is given by the 
condition 

a(0) = ku), (28) 

we can predict what will be its value at any moment. The parameter a(t) 
obeys the following Riccati differential equation 

a(t) = -ia 2 (t) + iu 2 . (29) 

Due to the Riccati procedure outlined before, we can linearize this equation 
by a substitution a = —in(t) /md(t) and then we know that the equations of 
the evolution of the parameters d(t) and n(t) are equivalent to the classical 
equations of motion (J8j) 

n(t) 

d(t) = h(t) = -muj 2 d(t). (30) 

m 
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Let us choose the initial conditions for n(t) and d(t) which are in accordance 
with (1281) as follows 



d(0) = -i\ — , n(0) = kVhmco. (31) 
V mu> 

From the solution of the classical equations of motion (|24|) one finds the 
time dependence of d(t) and n(t) 



d(t) = [cos(ut) — ik sinhdt)] , (32a) 

V muj 

n(t) = Vhmuj [kcos(iot) — isin(ut)] . (32b) 
It means that the solution of the Schrodinger equation has a form 

q>(x t) = ( * (^ COS ^ - iksin ^ \ 3 e »0(t) ™~ ( mujk cos M) - i sva.(ut) ^ 2 



\muj J \ kcos(ujt) — ism(u>t) J ^ \ 2H cos(tot) — iksm(ujt) 

(33) 

Such a wave function describes pulsating states known from the quantum 
mechanics courses. 

One should notice here that if we put k = 1 at the beginning of motion, 
then the time dependence will appear only in the phase of the wave function. 
It means that such a function is a stationary Gaussian wave packet of our 
system described by the Schrodinger equation (fl6|) . 



5. Gaussian stationary quantum states 

As we have shown above in the one dimensional example using the Ric- 
cati method we can not only predict the evolution of the shape but we 
can also find a stationary state of our system whose shape is constant in 
time. We will show now how to obtain such a stationary Gaussian state in 
a general N dimensional case. 



5.1. Conditions of stationarity 

Since the evolution equations of the parameters R(t) and P(t) are com- 
pletely separated from the rest, the Gaussian state can be stationary only 
when they are identically equal to zero. The condition that the shape is 
constant in time has the form 

! A - = °- < 34 > 

This condition guarantees that the normalization factor N(t) is constant in 
time, because it depends only on the shape of the wave packet. Therefore, 
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as it is seen from the equation (119dj) . the imaginary part of the matrix K is 
traceless. 



5.2. Solving the algebraic matrix Riccati equation 
To find a stationary Gaussian state 



exp 



— r ■ i(n • r 

2h 



we have to solve the following algebraic matrix Riccati equation 



= -iKl + iV 



(35) 



(36) 



where the matrix Kq describes the shape of a stationary Gaussian state. 

It is worth to notice that if we find two matrices D{t) and N(t) which 
satisfy the equations (|21|) and have a following form 



D{t) 
N{t) 



D -E(t), 
N -E(t). 



(37) 
(38) 



then the matrix Kq defined by the equation analogous to ([20 



K = --N ■ Dn 1 
m 



(39) 



is a solution of the equation ([36]) . 

The problem of finding such matrices is not a hard task. After all, we 
know that the columns of these matrices obey exactly the same equations 
as the classical positions and momenta. Therefore, if we take the classical 
eigenmodes as the columns of these matrices, then the matrix E{t) simply 
will be a diagonal matrix with the elements e %Wit and the matrices Do and 
N can be build from the amplitudes R L and Pj of the classical modes. 



5.3. Proper choice of classical eigenmodes 

The construction of the matrix Kq outlined above is not yet completed. 
Still it is not clear which N from 2N classical eigenmodes one should use 
in this construction. What are the requirements are necessary to make a 
good choice? The answer is: one should use such a set of classical modes 
that create a matrix Kq with a positive real part. It is a necessary re- 
quirement because we want this matrix to describe a real Gaussian, square 
integrable wave packet. Of course it can happen, that this construction is 
impossible, for example when one pair of characteristic frequencies is not 



10 



WAVE' FUNCTIONS PRINTED ON FEBRUARY 1, 2008 



real. In this case, our system does not have integrable quantum states but 
the construction of a stationary (unphysical) state is still possible. 

We should also notice here that even if the Hamiltonian is not bounded 
from below it can still be possible to construct integrable quantum station- 
ary states of our system. 



6. Other stationary states 

6.1. Evolution of wave packets with constant shapes 

It follows from the equations ([19]) that the motion of a Gaussian state 
with a constant shape K(t) = Kq, centered on the classical trajectory, is 
obviously possible. The wave function which describes such a situation has 
the form 



*(r,t) = Ne**®/*, 



exp 



m, „, - , n/ x . iP(t)-r 
- — {r - R(t))-K .(r - R(t)) + — ^- 



(40) 



where the vectors R(t) and P(t) obey the classical equations of particle 
motion (jHJ). Therefore, let us choose as a solution the most general physical 
trajectory generated by all possible modes of the system (fTZI) 



N 

R(t) = (AiBie*^* + A?fl i *e _iWi *) , (41a) 

i=l 
N 

P(t) = fcPi**** + Ai-P'e-***) . (41b) 
»=i 

Set of coefficients A, defines exactly one physical trajectory along which our 
gaussian solution with stationary shape moves. Wave function (|40p is just 
a multidimensional generalization of well known Glauber or coherent state. 

Without loss of generality one can assume that all modes used during 
construction of matrix Kq are labeled as R{ and Pj. All other modes are 
complex conjuncted to them. 

Because Ri and Pi are vectors used during construction matrix Kq 
therefore from (|39|) the following relations hold 

K -Ri = -—Pi, (42a) 
m 

K^-R* = i-Pj. (42b) 
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In addition from the equations (jlOj) follows that for the amplitude vectors 
Ri and Pi of any mode of the system following relations hold 



R* 




= Rj • P*, for 


i 3 


(43a) 


R* 


Pi 


= -Ri Pt 




(43b) 




Pi 


= m (^Cl + iu^j ■ Ri, 




(43c) 


im (uji + ujj) R{ 


Pj 


= Pi • Pj - m 2 Ri ■ V 


Rj, 


(43d) 


im (tOi — u>j ) Ri 


p* 

j 


= Pi ■ P* - m 2 R l ■ V 


T7* 

' j ' 


(43e) 



Using relations (|42h and (143)) one can easily show that the wave function 
Oj) can be represented in the following way 



tf(r,t) = Ne l 



where 



J e u exp 



N N 



w 



1 



N 



^ 6 X *i X *3 Ri-A-R^e-^^ 1 + K 

i=l j=l i=l 



r-A-R*e~ iuJlt 



*o(»0, 



A 

n 



K + K5, 



(44) 

(45a) 
(45b) 



and (fio is an arbitrary global phase factor. 

From the definitions of modes vectors Ri and R* it follows that any N 
from 2N such vectors is a base in the N-dimensional space. Therefore we 
can represent any N dimensional vector as a linear combination of chosen 
vectors. It means that there exist such coefficients Xi that 



r = J2xi(r)R*. 



(46) 

(47a) 
(47b) 

From the definition (I47ah one has the symmetry property jij = 7^. 
Using (|47|) one can represent the wave function of the coherent state 
as follows 



i=l 

To make next steps more clear let us introduce new quantities 

m 

^ = T 



A-R*, 



*(r,t) = iVe 4, 



J e _in °'exp 



N N N N 



i=l J=l 



i=l j=l 



*o(r) 



(48) 
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This form of the wave function allows us to find all stationary states of 
the system. 

6.2. Expansion into stationary states 

To get an expansion of the wave function (j48]) into stationary states 
(eigenstates of the Hamiltonian of the system) we should use the concept 
of, so called multidimensional Hermite polynomials. Properties of these 
polynomials are well studied [12] and a general information about them 
can be found in the Appendix [Bj The generating function of N dimensional 
Hermite polynomials lQ lt ...,n N (xi, ■ ■ ■ ,%n) has a form 



exp 



N N ^ 
i=l 3=1 



(ni,...,n N )=0 

(49) 

One can use this expansion in (|48p and then our wave function becomes 



(ni,...,n N )=0 



' N 

n 



_ e —irii(Jit 



H n 1 ,...,n N (zi>---^N) * (r). 



li=\ 

(50) 

Therefore our wave function is a superposition of the stationary states 
- the states whose evolution in time appears only in the evolution of the 
phase. For each set (n±, . . . ,n^) we have defined exactly one (up to the 
normalization factor) stationary state — the state which arise form the 
ground state by the excitations in the modes (i-th mode is excited rii times) 

*rn,...,n N (r) ~ K lt ...,n N ■ ■ ■ , x N ) *o(r). (51) 

It is worth to notice that the only place where the scale factors Aj appear 
in (|50p is the coefficient standing before the Hermite polynomial. Therefore, 
the wave function ([51]) does not depend on it. It means that set of the 
stationary states is independent on trajectory used during construction. 

It is obvious from the Schrodinger equation that the eigenvalue of the 

Hamiltionian in the state described by wave function (|5ip is equal H (Ya = i niuji + f^o 

All states with all possible sets of (n±, . . . , ton) form a base in the wave- 
function space. Because of the properties (|B.5P and (|B.7p of the multidi- 
mensional Hermite polynomials the completeness of this set is guaranteed 
and any quantum state of the system can be represented as a superposition 
of these states 



oo 

E 

(ni,...,n N )=0 



^(r,t)= V C nii ...,„ N (t)¥ nij .„ >nN (r). (52) 
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At this moment it is worth while to notice that the construction out- 
lined above uses effectively only this modes which were not used during the 
construction of the matrix Kq. One can easily see that amplitudes Ri do 
not play any role in the construction. It means that we have some kind 
of division of roles. Half of the modes is responsible for the creation of a 
positive defined matrix Kq and the ground state of the system. While the 
other half gives the stationary state of the system. 

This observation has a very nice physical interpretation. When the 
Hamiltonian of the system is bounded from below, one can show that the 
proper classical modes (which guarantee that the matrix Kq is positive de- 
fined) have positive frequencies. Therefore, when we create excited states, 
we have to use the modes with negative frequencies and their energy (eigen- 
value of the Hamiltonian) is growing with n^. On the other hand, when the 
Hamiltonian is not bounded from below, then some of the modes contained 
in the matrix Kq have negative frequencies. Therefore, when exciting one 
of these modes, one uses a mode with positive frequency and the energy of 
the generated states diminishes — it is not bounded from below. 

7. Conclusions 

In this paper we have shown how the classical-quantum correspondence 
for the linear systems is represented in the language of wave functions. The 
dynamics of Gaussian states is completely determined by the dynamics of a 
classical particle in a very unexpected way — the shape of Gaussian state 
evolves like a quotient of the classical momentum and position. This obser- 
vation enables us to construct a stationary Gaussian state of the system. In 
addition the classical dynamics of the center of the Gaussian state, leads to 
a complete set of stationary states. 
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Two mathematical theorems 
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Appendix A. 1 About the matrix Riccati equation 
The matrix Riccati equation in n dimensions has the form 
dk(t) 



dt 



aK\t) + K{t)-W + W L -K(t) + U. (A.l) 



where W and U are n x n matrices and a is a constant parameter. The 
Riccati theorem [5] says that K(t) is a solution of the equation fjA.lj) if 
and only if there exist matrices A(t) and B(t) which satisfy the following 
equations 

dA(t) ^ (j.^ t ) + w T -A(t), (A.2a) 



dt 
dB(t) 



aA{t) -W-B{t). (A.2b) 



dt 

and then the solution has the form 

K(t) = A(t)-B-\t). (A.3) 

It is obvious, that if there exist matrices A and B which satisfy the equations 
(|A.2[) then the matrix Kit) defined by an equation (|A.3P is a solution of 
the Riccati equation (|A.1|) . On the other hand, if K(t) obeys the equation 
([A~l]) then one can define the matrices A(t) and B(t) 



B(t) = Texp 



o 



aK(T) + W)dT 



(A.4a) 



A(t) = K{t)-B{t). (A.4b) 
One can then show that such matrices obey the equations (|A.2|) . 

Appendix A.2About the eigenvalues of some matrix 
Let us define 2n x In matrix M as follows 

M - (f I) (A.5, 

where W, A, B are real n x n matrices and A and B are symmetric. Then 
if A is an eigenvalue of the matrix (|A.5|) then —A is an eigenvalue of this 
matrix too. It means that the characteristic polynomial is of the order n 
in A 2 . This theorem follows from the observation that the matrices A4 and 
—M have the same set of eigenvalues because one has following reflexive 
property 

-M T =S.M-S-\ S=(° i f). (A.6) 
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Appendix B 

Multidimensional Hermite polynomials 

Multidimensional Hermite polynomials used in section 16.21 were intro- 
duced in [TT] and deep analysis of their properties was studied in [12]. Nowa- 
days these polynomials are often used in quantum optics contects [141 [15] 
in description of mixed and entangled states of electromagnetic filed. 

Multidimensional Hermite polynomials in N dimensions are functions 
of N variables x%,... ,xn- Let us assume that 7 is symmetric and positive 
N x N matrix. Then we define N dimensional Hermite polynomial of the 
order (m, ri2, ...,%) as follows 

(B.l) 

In one dimensional case it is a set of well known Hermite polynomials defined 
as 

H„(z) = (-l)^ 2 £U-* 2 . (B.2) 

In analogy to standard Hermite polynomials multidimensional Hermite 
polynomials can be also constructed with using a generating function which 
has a following form 



exp 



N N 

i=i 3=1 



m n N 
«1 



(ni,...,n N )=0 

(B.3) 

It is just a generalization of well known generating function of Hermite 
polynomials 



exp 



n=0 

Properties of the multidimensional Hermite polynomials are well studied 
[T2] . It was proved that they are orthogonal in the following sense 

N 

/ d N S e-^^H 7 lv .. inN (^)H7„ li ... imN ( S )=vr N / 2 n2^n,!^ mi . (B.5) 
Jul™ 

and any function f(x\, . . . of the class L"^, that is to say function for 
which the integral 

I d N xe- 1 2^\f(x 1 ,...,x„)\ (B.6) 
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is convergent, can be expanded in series of N dimensional Hermite polyno- 
mials |13j 

oo 

/(Xi, . . . , Xn) = fni,...,n N Hl u ...,n N (xi,- ■ ■ (B.7) 

(ni,...,n N )=0 

where 

/m,...,n N = ^ ^ J RN d * x e ~^' X f( x ^ • • • ^n) , nN (xi, . . . , x N ) . 

(B.8) 

Many other properties of the multidimensional Hermite polynomials may 
be found (cf. for example |12j). 
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